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A rich variety of Fermi systems condense by form- 
ing bound pairs, including high temperature 1 and 
heavy fermion 2 superconductors, Sr2Ru04 3 , cold atomic 
gases 4 , and superfluid 3 He 5 . Some of these form ex- 
otic quantum states having non-zero orbital angular mo- 
mentum. We have discovered, in the case of 3 He, that 
anisotropic disorder, engineered from highly porous silica 
aerogel, stabilizes a chiral superfluid state that otherwise 
would not exist. Additionally, we find that the chiral 
axis of this state can be uniquely oriented with the ap- 
plication of a magnetic field perpendicular to the aerogel 
anisotropy axis. At sufficiently low temperature we ob- 
serve a sharp transition from a uniformly oriented chiral 
state to a disordered structure consistent with locally or- 
dered domains, contrary to expectations for a superfluid 
glass phase 6 . 

Superconducting states with non-zero orbital angular 
momentum, L ^ 0, are characterized by a competi- 
tive, but essential, relationship with magnetism, strong 
normal-state anisotropy, or both 1 3 ' 5 . Moreover, these 
states are strongly suppressed by disorder, an important 
consideration for applications 7 and a signature of their 
unconventional behavior 3 ' 8 ' 9 . Although liquid 3 He in its 
normal phase is perfectly isotropic, it becomes a p-wave 
superfluid at low temperatures with non-zero orbital and 
spin angular momenta, L = S = l 10 . One of its two 
superfluid phases in zero magnetic field is anisotropic 
having chiral symmetry, where the handedness results 
from the orbital motion of the bound 3 He pairs about 
an axis t. This chiral superfluid, called the A-phase or 
axial state, is stable at high pressure near the normal-to- 
superfluid transition, Fig. la-c, while the majority of the 
phase diagram is the non-chiral B-phase, having isotropic 
physical properties. The stability of the A-phase is at- 
tributed to strong-coupling effects arising from collisions 
between 3 He quasiparticles 10 . However, in the presence 
of isotropic disorder these strong-coupling effects are re- 
duced and the stable chiral phase disappears 11 ' 12 , Fig. la. 
Here we show that anisotropic disorder can reverse this 
process and stabilize an anisotropic phase over the entire 
phase diagram, Fig. lc. 

For many years it was thought to be impossible to 
introduce disorder into liquid 3 He since it is intrinsi- 
cally chemically and isotopically pure at low tempera- 
tures. Then it was discovered 15 ' 16 that 3 He imbibed in 
~ 98% porosity silica aerogel, Fig. Id, is a superfluid 
with a sharply defined transition temperature 12 , but re- 
duced from that of pure 3 He. To test predictions that 
isotropic disorder favors isotropic states and conversely 17 , 
we have grown a 97.5% porosity anisotropic aerogel with 
growth-induced radial compression 13 , effectively stretch- 



ing it along its cylinder axis by 14.3%. Experiments using 
uncharacterized stretched aerogels have been previously 
reported 18 ' 19 and are in disagreement with the work pre- 
sented here. Silica aerogels, as in Fig. Id, are formed 
by silica particles « 3 nm in diameter, precipitated from 
a tetramethylorthosilicate solution, and aggregated in a 
diffusion limited process. After supercritical drying we 
obtain a cylinder as shown in Fig. Id. Our numerical 
simulation of this process indicates that the gel particles 
form strands, see Fig. Id. According to theory 20 , t is con- 
strained to point outward from the strand. Stretching an 
aerogel tends to align the strands, thereby forcing the chi- 
ral axis, on average, to be in an easy-plane perpendicular 
to the aerogel anisotropy axis 6 . On the other hand, for 
an axially compressed aerogel t should be aligned with 
the compression axis 6 . 

We use pulsed nuclear magnetic resonance (NMR) 
(Methods Section), to identify the 3 He superfluid state 
and determine the magnitude of its order parameter, A. 
The equilibrium nuclear magnetization is tipped by a 
pulsed magnetic field to an angle f3 away from the static 
applied magnetic field, H, and the Fourier transform of 
its free precession is the NMR spectrum. The frequency 
shift of the spectrum, Acj, from the normal- state value 
centered at the Larmor frequency, ujl , is a direct measure 
of the temperature dependent order parameter, A(T). 
For the A- and B-phases of 3 He 10 the frequency shifts 
are, 

A-phase: 

a -<«=H(H c «<») 

B-phase: 

Auj b (P) = /3<104° 

0| = A|( 6 ^^), (2) 

where Vt\ B is the square of the longitudinal resonance 
frequency, which is proportional to A\ B and determined 
by constants of the normal fluid: gyromagnetic ratio, 7; 
dimensionless dipole coupling constant, A^; single spin 
density of states at the Fermi level, N F ; as well as the 
magnetic susceptibilities xa and xb- The A-phase sus- 
ceptibility is temperature independent and equal to the 
normal-state value, xn', however, xb is temperature de- 
pendent and suppressed relative to xn- 
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FIG. 1. Superfluid phase diagrams, a-c) Comparison of phase diagrams of 3 He in silica aerogels with various types of 
anisotropy in the limit H — ► 0. High spatial uniformity of the aerogels, better than 99.97%, was established using optical cross 
polarization methods 13 for a) and c). The solid black curves are fits to T ca (P) using the GL-scattering theory (Supplementary 
Information) and the solid (dashed) red curves correspond to the superfluid transition T c (B — >• A transition Tba) for pure 
3 He at H — 0. For subscripts, the lowercase a signifies aerogel and uppercase A and B represent the corresponding phases of 
superfluid 3 He. a) A uniformly isotropic 98.2% porosity aerogel 12 , b) A 98.0% porosity aerogel axially compressed by 10% from 
torsional oscillator experiments by Bennett et al. 14 . In both a) and b) the B-phase occupies the majority of the P — T plane, 
c) For the 97.5% porosity stretched aerogel used in the present work the A-phase is stable throughout the phase diagram. The 
open squares (Td) mark the transition to a disordered chiral phase and the open circle at T — indicates a quantum critical 
point at which the superfluid order parameter vanishes (see Fig. 2b). d) Photograph of an aerogel with ~ 98% porosity, 1 cm 
in diameter, partially extracted from the glass tube in which it was grown. The microstructure of aerogel based on numerical 
simulation of diffusion limited aggregation of silica particles, with diameter 5, into a network of strands with most probable 
spacing, £ a . The blow-up sketch shows the favorable directions (red arrows) of the chiral axes, £, perpendicular to the strand. 



For 3 He in our anisotropic aerogel we observe a sharp 
transition to a superfluid state at a temperature, T ca , 
marked by the onset of NMR frequency shifts for small 
tip angle, = 8°, Fig. 2a. Our comparison of the data 
with either of the two pure 3 He states based on Eq. 1 
and 2, suggests that the superfluid is an axial p-wave 
state like the A-phase. In the Ginzburg-Landau (GL) 
limit, T < T c , the initial slope of Auja is proportional to 
the square of the order parameter magnitude A\ where 
A 2 A (T) = A 2 A0 (1 -T/T c ). We find in our anisotropic 
aerogel that A\ has a linear pressure dependence, as 
shown in Fig. 2b, just as is the case for pure 3 He-A 21 ' 22 . 
It is, however, substantially reduced in magnitude and 
extrapolates to a critical pressure of P c = 7.9 bar at 
T = where the order parameter vanishes at a quan- 
tum critical point 23 . We have analyzed the transition 
temperatures T ca (P) using GL theory 24 (Supplementary 
Information) to determine the mean free path A = 113 
nm and the silica particle correlation length £ a = 39 nm. 
Our calculated phase diagram compares very well with 



the data in Fig. lc. With the same parameters we have 
also calculated A^q (Supplementary Information), shown 
in Fig. 2b, and find excellent agreement with the order 
parameter determined from the frequency shifts, Eq. 1, 
consistent with a suppressed A-phase. Moreover, extrap- 
olation of Auj(P) is at a pressure that coincides 
with the critical pressure from the calculated phase dia- 
gram, demonstrating important consistency between ex- 
periment and theory. 

The tip angle dependence of the NMR frequency shift 
is a fingerprint of a specific p-wave state 10 . Our measure- 
ments of Auj((3), Fig. 2c, follow the expected behavior of 
the chiral axial state 26 , Eq. 1. The magnitude of the shift 
is reduced compared to that of pure 3 He- A as shown in 
Fig. 2b. Furthermore, the axial state belongs to the class 
of equal spin pairing (ESP) states, having the same sus- 
ceptibility as the normal Fermi liquid. We have directly 
demonstrated that this is the case as shown in Fig. 2d. 

Lastly, we have looked for a possible phase transition to 
some other p-wave state by cooling to our lowest temper- 
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FIG. 2. Identification of the superfluid state, a) The frequency shift of the NMR spectrum for a small tip angle, /3 = 8°, 
as determined from gaussian fits. The superfluid transition temperature at different pressures is marked by the onset of positive 
frequency shifts. The NMR spectra have been corrected for fast exchange of liquid spins with several layers of paramagnetic 
solid 3 He on the aerogel surface 16 ' 25 , b) The square of the order parameter magnitude as a function of pressure in our aerogel 
(closed green circles) was calculated from the initial slopes in a) using Eq. 1, showing the suppression relative to pure 3 He- A 
(open and closed squares) 21 ' 22 . The experimental results are compared to the predictions of the GL-theory (open blue circles) 
(Supplementary Information). The critical pressure, P c — 7.9 bar is determined from an extrapolation (dashed line) of A\ 
to zero, indicated by the arrow. The solid line is a linear fit to the pure A-phase data, c) Identification of the superfluid as 
the axial p-wave state follows from the agreement with theory 26 (dashed blue curve) for the dependence of Auj on tip angle, 
Td < T < T c . d) The susceptibility for stretched and isotropic aerogel samples taken on warming at a pressure P c± 26 bar 
in H = 196 mT after subtraction of a paramagnetic background from solid 3 He 16 . For the isotropic aerogel the field-induced 
A-phase region has a constant susceptibility below T ca and a very sharp jump at the BA-transition 12 . A similar jump in 
susceptibility would be anticipated for the stretched aerogel if the 3 He were in the B-phase or any other non-ESP phase. 



atures, < 650 /iK, in a substantial magnetic field of 196 
mT, and measuring x an d Aco> on warming. For exam- 
ple, if there had been a transition to a B-phase on cooling 
we would have observed a discontinuous increase of \ on 
warming from the B-phase to a magnetic field-induced A- 
phase, as is characteristic of the first order B-to-A transi- 
tion in an isotropic aerogel (Fig. 2d) 12 . Correspondingly, 
we would also have observed a discontinuous change in 
Auj. The absence of these discontinuities demonstrates 
that an ESP superfluid state in stretched aerogel is stable 
down to our lowest temperatures. 

Our three independent NMR measurements show that 
anisotropic stretched aerogel stabilizes anisotropic super- 
fluidity throughout the phase diagram, Fig. lc. This 
is in stark contrast to 3 He in the presence of isotropic 
disorder 11 ' 12 , Fig. la, or for anisotropic disorder estab- 
lished by axial compression, Fig. lb 14 . In both of these 



latter cases the B-phase is dominant. Understanding the 
stability of chiral superfluidity, stabilized by anisotropy, 
is an interesting open problem, in part because topolog- 
ically non-trivial pairing states have unusual structure. 

Our NMR results reported in Fig. 2a-c were restricted 
to temperatures within ~ 20% of T ca . In this high tem- 
perature region there is negligible increase in the NMR 
linewidth as compared to the normal fluid, see Fig. 3, in- 
dicating that the chiral axis, i, has a well-defined direc- 
tion orthogonal to both the aerogel anisotropy axis 6 and 
the magnetic field. This orientation minimizes the mag- 
netic dipole energy producing the so-called dipole-locked 
condition on which the validity of Eq. 1 depends 10 , and 
with which our tip angle results agree. The combined 
effects of the dipole energy and aerogel anisotropy define 
a uniform direction of £ for the A-phase of superfluid 3 He 
in our stretched aerogel. 
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FIG. 3. The disorder transition. Below the disorder transition, T < Td (vertical dashed lines), the NMR frequency shift 
drops and the linewidth abruptly increases. The frequency shift is defined as the first moment of the NMR spectrum and the 
linewidth is the full width at half maximum. Both quantities have been corrected for the effect of fast exchange between liquid 
and surface solid 3 He spins 16 ' 25 . Data were taken on warming. 



At lower temperatures, we have discovered a different 
phase characterized by a significantly broadened NMR 
line. The data indicate that this superfluid state is no 
longer the homogeneous A-phase with £ _L H. On warm- 
ing there is an abrupt increase in the first moment of 
the NMR spectrum (frequency shift) and a correspond- 
ing decrease in the linewidth at a well-defined temper- 
ature, Td (Fig. 3). The transition at Td is independent 
of magnetic field for 31.1 < H < 196 mT. In addition, 
the magnetic susceptibility does not change at this tem- 
perature; however, Td is pressure dependent as shown by 
the open squares in the phase diagram of Fig. lc. We 
interpret these results as evidence for a disordered ESP 
phase with a possible phase transition at Td. 

The frequency shift is a measure of A^, given in Eq. 1 
for the dipole-locked case where the angle between the 
chiral axis and the magnetic field is = cos _1 (^ • H) = 
90°. However, in general the shift also depends on 6 
which for small tip angles Au(0) reduces to 27 , 



Q2 

Acj(0) = ^cos(20). 

2ujl 



(3) 



This allows us to extract the angular distribution, P(#), 
of £ relative to H, from the NMR spectra both above and 
below Td. The observed spectrum is given by the con- 
volution product of the normal-state lineshape and the 
frequency distribution P(oj) that corresponds to P(0) de- 
termined by Eq. 3 (Supplementary Information). Above 
Td the NMR spectrum (Fig. 4a) can be fit with an £- 
distribution (Fig. 4b) that has a single component cen- 
tered at = 90°, as sketched in Fig. 4c, with angular 



spread AO = 13° ± 2.3°. Examination of the spectrum 
below Td, Fig. 4a, indicates that it is bimodal, suggest- 
ing that P(9) has three main components, which we take 
to be gaussian functions with adjustable position, width, 
and relative weight. From our fits we determine these 
parameters and obtain P(0), Fig. 4b. For the spectra 
below Xd, the disordered state is composed of three do- 
mains. Approximately 1/3 of the sample remains with 
£ JL H while the majority, ~ 2/3 of the distribution, 
has nearly zero frequency shift with = 44 ± 0.5°. The 
theory, Eq. 3, does not distinguish between angles that 
are symmetric about 6 = 90°, and thus we assume that 
P(6) = P(180° - 0) as shown in Fig. 4d. Our inter- 
pretation of the disordered state in terms of distinct £- 
domains is different from the two-dimensional "orbital 
glass" phase defined by a random distribution of £ con- 
fined to a plane perpendicular to the stretching axis 6 . It 
is noteworthy that the possible existence of orbital do- 
mains, such as we observe here, has been central to the 
discussion of time-reversal symmetry breaking in the un- 
conventional superconducting state in S^RuC^ 28 . 
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FIG. 4. Distribution of the direction of the chiral axis, a) The convolution product of the normal-state NMR line (red 
curve) with a best-fit angular distribution of the chiral axis, P(0), gives an excellent representation (dashed black curves) of 
the NMR spectra for the superfluid state above (left panel, bold green curve) and below (right panel, bold green curve) the 
disorder transition (Supplementary Information). The spectra were obtained with H — 31.1 mT. b) The angular distributions 
for the chiral axis, £, determined from a). The dashed pink curve is the distribution for Td < T < T ca . c) Stretching the aerogel 
along the cylinder axis forces t into an easy-plane perpendicular to the direction of strain 6 . To minimize the dipole energy t 
must be perpendicular to H, as illustrated, d) According to our analysis the distribution splits into three ^-domains for T < Td- 



I. METHODS 

In our experiments the external magnetic field H was 
31.1, 95.5, or 196 mT and oriented perpendicular to 
the aerogel cylinder axis. The sample was 4.93 mm 
long and had a diameter of 3.43 mm. The susceptibil- 
ity was obtained by numerically integrating the phase- 
corrected absorption spectrum, while the frequency shift 
and linewidth were determined from the power spectrum. 
The sample was cooled using adiabatic nuclear demag- 
netization of PrNis. NMR on 195 Pt was used for ther- 
mometry at H = 95.5 and 196 mT and calibrated rela- 
tive to the known phase diagram of pure superfluid 3 He 



from a volume of liquid outside of the aerogel sample 
equaling 30% of the total liquid. For data acquired at 
H = 31.1 mT, temperature was determined from the 
Curie- Weiss dependence of the solid 3 He adsorbed to the 
aerogel surface 16 ' 25 . 
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SI.1. GINZBURG-LANDAU THEORY 

The Ginzburg-Landau (GL) free-energy functional for £>-wave, spin-triplet pairing of 3 He 
quasiparticles in a homogeneous isotropic medium is defined in Ref. [SIl] in terms of the 
3x3 matrix order parameter, A$i, which transforms as a vector under spin rotations with 
respect to index 5, separately as a vector under orbital rotations with respect to index z, 



The material coefficient, a(T) ~ a'{T — T ca ), depends on the single spin density of states 
at the Fermi surface, N F) the pairing interaction, and the effects of scattering by the ran- 
dom impurity potential, all of which determine the transition temperature, T ca . The prime 
denotes differentiation with respect to temperature. The fourth-order material parameters, 
= 1...5}, which are also sensitive to impurity scattering, determine condensation 
energy and the structure and residual symmetry of the superfluid condensate. 

The homogeneous isotropic scattering model (HISM) of Thuneberg et al [SIl] was ex- 
tended by Sauls and Sharma in Ref. [SI2] to incorporate the effect of correlations between 
silica strands and clusters. In this model, the aerogel correlation length, £ aj limits the pair- 
ing transition for small superfluid coherence lengths, ^ j p/27rfcsT ca < £ a (high pressures), 
while the mean free path, A, determines the transition when the superfluid coherence length 
is larger than the aerogel correlation length. This condition is achieved at low pressures, 
and always in the limit T ca — > 0. Indeed this limit can be used to determine A from the 
critical pressure at which T ca (P c ) = 0. The two pair-breaking length scales define an ef- 
fective pair-breaking parameter, x — x/(l + Ca/ X )> (a = £a/A where x = hvp /^ksTX is 
the depairing parameter in the HISM [SI3] and vp is the pressure dependent Fermi velocity 
of 3 He [SI4]. The superfluid transition temperature in aerogel, T ca , is determined by the 
condition a(T ca ) = in the following equation, 



where T c is the transition temperature of pure 3 He. In the limit x = 0, Eq. SIl reduces to 
the case for pure 3 He, namely a = (l/3)iVp ln(T/T c ). 
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FIG. SI1. Determination of Scattering Parameters Aerogel correlation length (red circles) 
and standard deviation of £ a as a function of A (blue squares). The minimum in the standard 
deviation corresponds to A = 113 nm and ^ a = 39 nm. 

The correlation length is a pressure independent length scale of the aerogel microstruc- 
ture. We have used Eq. 2 to fit our measurements of T ca (P) in the stretched aerogel (Fig. lc) 
with the constraint that ^ a be pressure independent. From this fit we find A = 113 nm and 
£ a = 39 nm. In Fig. SI1 we present the standard deviation of t; a (P) as a function of A. 
Using these same values of A and £ a we have calculated the square of the axial state order 
parameter amplitude in aerogel from the GL theory, 



A 2 



^ 2 a' {Tea) (3 A 



(SI3) 



M T c a'(T c ) /3a 

In Eq. SI2 A 2 A0 is the value for the pure A-phase [SI5, SI6], and /3a (/3a) is the combination of 
/5-parameters appropriate to pure (aerogel) 3 He-A, i.e. (3a = /?2 + /?4 + /35 (Pa — ^2 + ^4 + ^)- 
The pure /5-parameters were taken from Ref. [SI7]. For calculating the aerogel ^-parameters, 
we have extended the results of Thuneberg et al. [SI1] for the HISM, to include the combined 
effects of the aerogel correlation length and mean-free path using the correlated scattering 
model of Ref. [SI2]. 

SI.2. NMR LINE SHAPE ANALYSIS 



We have performed an analysis of the NMR line shape to obtain the distribution of the 
chiral axis, £ : in our stretched sample. A distribution in the direction of £ relative to the 
external magnetic field, will lead to a distribution in the intrinsic superfluid dipole field 



and hence a distribution of frequency shifts in the NMR spectrum. If 9 = cos _1 (^ • H) is 
the angle between the chiral axis and the magnetic field, then for small tip angle, /5, the 
frequency shift in the A-phase reduces to [SI8] , 

Q2 

Auj a (9) = ^cos(2#). (SI4) 

For example, if the chiral axis orientation is unique and the normal state lineshape is given 
by F n (uj — ujl) then the lineshape in the superfluid, F S) is unchanged except for a constant 
shift, i.e. F s = F n (uj — uj l — Auja) where Acja is given by Eq. SI3. However, for an angular 
distribution of t, P{9\ there is a distribution of frequency shifts, P((jS), determined by 
Eq. SI3. This results in a more complex NMR lineshape that is given by the convolution 
product of the normal state line and P{(jo) ) 

F s (uj) = J P{J)F n a/) duo'. (SI5) 

Since the normal state and superfluid spectra are obtained experimentally, P{6) can be 
determined by fitting the spectrum in the superfluid phase with Eq. SI4. 

In general an inversion problem, like this one, requires constraints which make it difficult 
to determine the uniqueness of the distribution. We choose the model that the ^-distribution 
can be represented by a sum of gaussian functions with adjustable position, i? width, 
and relative weight, 

P{9) = Y,Aie^)\ (SI6) 

i 

and fP (6) d6 — 1. In Fig. 4a we present NMR spectra (bold green curves) and correspond- 
ing fits (dashed black curves) for temperatures above (left panel) and below (right panel) 
the disorder transition Td with the results for P{6) shown in Fig. 4b. For Td < T < T ca 
we find that the spectrum is well-represented with a distribution centered at 9 = 90°, cor- 
responding to the dipole-locked configuration, with a spread of 13.0° ± 2.3°. For T < Td, 
we find that the ^-distribution is composed of domains as described in the text. Approxi- 
mately 1/3 of the sample remains with t perpendicular to H with a spread < 4° as shown 
in Fig. SI2. where we present the squared residual of the fit as a function of the width of 
the dipole-locked peak while keeping the other parameters fixed. The remaining ~ 2/3 of 
the superfluid has 9 = 44° ± 0.5° with a width of 4.9° ± 0.6° as shown in Fig. SI3. This 
distribution fits our spectrum well for T < T d . To check its uniqueness we have attempted 
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FIG. SI2. Squared residual versus the width of the dipole-locked peak for T < 
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FIG. SI3. Squared residual versus the a) position, b) percent weight, and c) width of the dipole- 
unlocked peak for T < 



to add an additional peak with variable relative weight, w v , at the position 9 = 23°, with 
a width of 9 = 6° and we find that the best fit to the spectrum corresponds to w v = 
as shown in Fig. SI4a. To emphasize this point, in Fig. SI4b we show the lineshape that 
results for the case w v = 15% that clearly adds substantial weight in the spectrum which 
is inconsistent with our experiment. We infer, therefore, that the disordered state has an 
^-distribution which is non-monotonic with three principal components as shown in Fig. 4b. 
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FIG. SI4. Attempt to add an additional peak a) Squared residual versus the weight of an 
additional peak at 9 — 23°. b) Fit (dashed black curve) to the NMR spectrum below (bold 
green curve) produced with w v = 15%. The normal state line is presented as the fine red curve 
and the fit is performed using Eq. SI4. Significant discrepancy between the best fit and the NMR 
spectrum indicates that extra weight in the ^-distribution between 44° and 90° can be ruled out. 
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